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Abstract 

We compute the average polarisation asymmetry from the Klein-Nishina differ¬ 
ential cross section on free electrons at rest. As expected from the expression 
for the asymmetry, the average asymmetry is found to decrease like the inverse 
of the incident photon energy asymptotically at high energy. We then compute 
a simple estimator of the polarisation fraction that makes optimal use of all the 
kinematic information present in an event final state, by the use of “moments” 
method, and we compare its statistical power to that of a simple fit of the az¬ 
imuthal distribution. In contrast to polarimetry with pair creation, for which 
we obtained an improvement by a factor of larger than two in a previous work, 
here for Compton scattering the improvement is only of 10-20 %. 
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1. Cosmic-source polarimetry: the high-energy frontier 

Polarimetry is a powerful diagnostic of specific phenomena at work in cos¬ 
mic sources in the radio-wave and optical energy bands, but very few results 
are available at high photon energies: the only significant observation in the 
X-gamma energy range, to date, is the measurement of a linear polarisation 
fraction of P = 19 ± 1% of the 2.6 keV emission of the Crab nebula by a 
Bragg polarimeter on board OSO-8 [T]. At higher energies, hard-X-ray and 
soft-gamma-ray telescopes that have flown to space in the past (COMPTEL 
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P], BATSEPj) were not optimized for polarimetry, and their sensitivity to po¬ 
larisation was poor. Presently active missions (Integral IBIS [4] (5 i and SPI j 6 [) 
have provided some improvement, with, in particular, mildly significant mea¬ 
surements of P = 28 ± 6 % (130 to 440 keV 0) and P = 47+] 3 % (200 to 800 
keV |5) for the Crab Nebula. A number of Compton polarimeter/telescope 
projects have been developed, some of which also propose to record photon con¬ 
versions to e + e _ pairs. A variety of technologies have been considered, such 
as scintillator arrays (POGO [7], GRAPE [8], POLAR [9]), Si or Ge microstrip 
detectors (MEGA [TO], ASTROGAM [IT]) or combinations of these (Si + LaBr 3 
for GRIPS H2], Si + CsI(Tl) for TIGRE [Til]), semiconductor pixel detectors 
(CIPHER [H|) and liquid xenon (LXeGRIT [T5]) time projection chambers 
(TPC). In most Compton telescopes the reconstruction of the direction of the 
incident photon provides an uncertainty area which has the shape of a thin cone 
arc. The tracking of the recoil electron from the first Compton interaction with 
a measurement of the direction of the recoil momentum, as is within reach with 
a gas TPC, allows to decrease the length of the arc and therefore to improve 
dramatically the sensitivity of the detector (HU and references therein). 

Some of these telescopes are sensitive to photon energies up to tens of MeV 
in the Compton mode, but their sensitivity to polarisation above a few MeV is 
either nonexistent or undocumented. 

2. Polarisation asymmetry and average polarisation asymmetry 

As is well known, the sensitivity to polarisation of Compton scattering is ex¬ 
cellent at low energies (Thomson scattering), as the polarisation asymmetry A, 
also known as the modulation factor and defined by the phase-space dependence 
of the differential cross section, 

^ a [1 + AP cos [2(0 - 0 O )]] , (1) 

reaches —1 at a polar angle 6 of 90° (Fig. 1 of Ref. [T7]). In this expression, 0 
is the azimuthal angle, that is the angle between the scattering plane and the 
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direction of polarisation of the incident photon. Unfortunately, A is decreasing 
with energy, and as the precision of the measurement scales as ap oc 1 /(A\fN) 
when the background noise is negligible and where N is the number of signal 
event, the sensitivity of Compton polarimetry decreases at high energies. With 
the goal of a quantitative assessment of this sensitivity, in this paper we compute 
the average polarisation asymmetry (*4) from the Klein-Nishina differential cross 
section on free electrons at rest |T 8 | [T5j. (A) is defined from the differential cross 
section in </>, that is after the full differential cross section (eq. ®) has been 
integrated over the other variables that describe the final state: 


da 

d 4> 


oc [1 + (A)P cos [2(cf> - </>o)]]. 


( 2 ) 


Following Heitler J2D], the doubly differential cross section for linear polarised 
radiation reads: 


da rg 2 
dU 2 


-2 sin 2 9 cos 2 < 


(3) 


where x = fc/fco, fco and k are the energy of the incident and scattered qs, 
respectively; 8 is the scattering angle, that is the polar angle of the direction of 
the scattered 7 with respect to the direction of the incident 7 . The differential 
element dS2 is sinf?d0d0 as usual. In the case of partially polarised emission 
with polarisation fraction P, the differential cross section becomes: 


da = '4x 2 


1 


x H-sin" 9 (cos (2 (j))P + 1) 

x 


sin 9<A9A(j>. 


(4) 


The minus sign reflects the fact that photons Compton scatter preferentially 
into the direction perpendicular to the orientation of the electric field of the 
incoming radiation. The energy of the scattered 7 is related to 8 from energy- 
momentum conservation: k = fco/[l + fc 0 (l — cos 0 )], x = 1/[1 + &o(l — 
cos#)], cos 9 =1— (1/x — 1) /fco, sindd# = — dx/(koX 2 ) and sin 2 9 = 
2 \\/(xko) — 1/fco] — \l/(xko) — 1/fco] 2 . We then obtain [20112T] : 
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k varies in a range such that —1 < cos0 < 1, that is 1/(1 + 2fco) < x < 1. 
The distributions of these kinematic variables are shown in Fig. [lj After an 
elementary integration over x, we obtain: 


da 


2 1 + ko f 2fco(l + ko) 

° I ^0 V l + 2A:o 


In (1 + 2ko 


ln(l + 2fc 0 ) 

2fc 0 


1 + 3ko (2fco — {ko + 1) log(2fco + 1)) 

- 0T2W +- tl - cos (2 « p 




( 6 ) 


that is a total cross section of mmm-- 


a = 27rrg x (7) 

l + ^o /2fco(l + fco) , f i , oll A ln(l + 2fco) 1 + 3fco 

[ntr V 1 + 24 " 1,1 (1 + 2ko) ) + 2t 0 - (TT^fJ ■ 


Equating the constant term and the term proportional to cos (2 cf>)P in eqs. 
([2| and ([d]), we obtain for the average polarisation asymmetry: 


(A) 


(2 /cq ~ (ko + 1) log(2fco + 1)) 


(1 + kg) 
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1 + 2/cq 
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( 8 ) 
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Figure 1: Spectra of the azimuthal and polar angles (p and 9 , of cos0, of the fraction x of the 
incident photon energy carried away by the scattered photon, and of the ID and 2D weights 
w and i^opt, for incident photon energies O.lmc 2 , me 2 , and 10mc 2 , all for a fully polarized 
beam. 

We now examine two limiting cases: 
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Figure 2: Absolute value of the average asymmetry in Compton scattering on free electrons 
at rest, as a function of the incident photon energy in electron rest-mass units. Thick solid 
line: full expression (eq. (|8|). Dashed line: high-energy approximation (eq. 0 ) . Thin solid 
line: high-energy asymptote. 


• At low energies, k 0 « 0, eq. ([ 6 ]) reduces to: 

do- = ^ [1 — cos (2c/)P/2] y dcj>, (9) 

which results in a total cross section of a = 8 rQ 7 r/ 3 , i.e., the Thomson 
cross section. The low-energy average asymmetry is (A) = —1/2. 


At high energies, 


do- = -f- 
2kn 


log( 2 fc 0 ) + l 


4 — 2 log(2/c 0 ) 

P-:- COS [Zip) 


d </, 


which results in a total cross section of a = 7 rro(log( 2 £;o) + \)/ko. 
high-energy average asymmetry is 


( 10 ) 

The 


(A) 


4(log 2 fc 0 - 2 ) 
fc 0 ( 21 og 2 fc 0 + 1 )' 


( 11 ) 


The average asymmetry decreases at high energies, asymptotically approach¬ 
ing ( A) ~ — 2/fco- The variation of the average polarisation asymmetry of pho¬ 
ton Compton scattering on free electrons at rest (eq. ([ 8 ])) is compared to its 
high-energy approximation (eq. @) in Fig. § The absence of sensitivity of 


5 







Compton polarimeters at high energies [22] is due to this strong decrease of 

m\. 


3. Optimal variable for polarisation measurement 


The value of the polarisation fraction P is classically obtained by a fit to 
the (j> distribution. A way to improve the polarisation sensitivity is to make an 
optimal use of the information contained in the multi-dimensional probability 
density function (pdf) through the use of an optimal variable ([23] and references 
therein), that is, of a weight w(S 1) such that the P dependence of the expectation 
value E(w) of w allows a measurement of P, and that the variance of such a 
measurement is minimal. The solution, up to a multiplicative factor, is (eg. 


my 

^opt 


<91np(fi) 

dP 


In the present case of a polarisation measurement: 


( 12 ) 


p{n) = f(n) + p x g(n), 


(13) 


with f /(fl)dfl = 1 and f g(fl)dfl = 0, we obtain: 

g(fl) 

° pt fm + Pxg(Q)’ 
that is, if |P x g(fl)| is small compared to 


^opt — 


g{fy 

/(«)’ 


(14) 


(15) 


P 


The 1 st moment of w opt is (w opt ) = j y x [/(fl) + P x g{ fl)] df2 = 

/• 9 2 m 


m) 


dfl, which is proportional to P. The expressions for / and g are 


obtained from the measured values of <f> and 9 (and therefore of x) by equating 
the constant term and the term proportional to cos (2 <j>)P in eqs. (13) and ([l]). 
The spectrum of w opt is shown in Fig. [l] for a fully polarised beam. We can 
see that |ui opt | is most often much smaller than unity (beware the vertical log 
scale), so that our neglecting |Pxg(f2)| in the expression of f(Si) + Pxg(Q) was 
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legitimate. The asymmetry, the non-evenness of the Wopt distribution makes the 
non-zero average due to the beam polarisation explicit. 

Moment’s methods are equivalent to a likelihood analysis in the case where 
the pdf is a linear function of the variables that one aims to measure, as is the 
case here, but they are much simpler to instantiate as one just has to compute 
w 0p t(6,(j),x), and average it over the whole statistics. Although the analysis of 
experimental data is beyond the scope of this paper, the following considerations 
apply: 

• Background subtraction reduces to a simple subtraction in computing the 
average of w op t . Their n-dimensional parametrization is not needed. 


• Likelihood methods need the use of a n-dimensional parametrization of the 
acceptance, or efficiency, for correction. This is pretty simple in the case of 
Compton scattering for which the final state is described by only two vari¬ 
ables, but for higher-dimensional systems, producing enough Monte Carlo 
(MC) statistics and determining a parametrization becomes a nightmare: 
in that case the use of a moments-based efficiency correction becomes 
mandatory (for a real-case presentation see eg., section IV.A, eqs. (18)- 
(24) and VLB eqs. (47)-(49) of [23]). 


In the “reduced” ID case of eq. Q, w opt becomes w = 2cos2</> and the 
estimator for (A)P is (w) [23] . The uncertainty then reads: 


(16) 

that is, in the case of Thomson scattering ((M) = —1/2), ap ~ a/(8 — P 2 )/N. 
Needless to say, in the case where the direction of the polarisation of the emis¬ 
sion of a particular cosmic source “on the sky” is unknown, a combined use of 
(2 cos 2<j>) and of (2 sin 2 <j>) should be used. 

The performance of the 2D estimator (w op t) is compared to that of the ID 
(w) by the comparison of the ratios of the RMS width normalized to the mean 
value: 

RMS,„ 0 pt /RMS W (1?) 


(^opt) 


”> ' 
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Figure 3: Ratio r of the figures of merit of the 2D to ID estimators of the linear polarisation 
fraction, as a function of the incident photon energy in electron rest-mass units. 


In contrast with polarimetry performed with e + e - telescopes, for which 
an improvement in the precision of the measurement of the linear polarisation 
fraction by a factor of larger than two is at hand (Fig. 21 right of [23]), in the 
case of Compton polarimeters the improvement is found to be marginal, varying 
from S3 20% at low energy to « 10% at high energy (Fig. |3|. These results are 
in qualitative agreement with those obtained at 100 keV by a likelihood analysis 
of the doubly differential cross section □ 21 - 

In summary, we have obtained the expression for the average polarisation 
asymmetry, or modulation factor, of Compton scattering on free electrons at 
rest, eq. ([8]), Fig. [2j We have then obtained a simple optimal estimator of the 
polarisation fraction P that makes use of all the information (azimuthal and 
polar angles of the scatter), avoiding the technicalities of a maximum likelihood 
analysis but with the same performance. 
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